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Abstract 

Icing scaling tests in the NASA Glenn Icing Research Tunnel (IRT) were performed on swept wing 
models using existing recommended scaling methods that were originally developed for straight wing. 
Some needed modifications (Ref. 1) on the stagnation-point local collection efficiency (i.e., p 0 ) calculation 
and the corresponding convective heat transfer coefficient for swept NACA 0012 airfoil models have been 
studied and reported in 2009, and the correlations will be used in the current study. The reference tests 
used a 91.4-cm chord, 152.4-cm span, adjustable sweep airfoil model of NACA 0012 profile at velocities 
of 100 and 150 knot and MVD of 44 and 93 pm. Scale-to-reference model size ratio was 1:2.4. All tests 
were conducted at 0° angle of attack (AoA) and 45° sweep angle. Ice shape comparison results were 
presented for stagnation-point freezing fractions in the range of 0.4 to 1 .0. Preliminary results showed that 
good scaling was achieved for the conditions test by using the modified scaling methods developed for 
swept wing icing. 


Introduction 

The results presented here are part of an effort to develop appropriate scaling methods for swept wing ice 
accretion. Previously, Anderson (Ref. 2) has completed a NASA report that gave a detailed technical review 
of recommended scaling methods for ice accretion on unprotected, unswept aerodynamic surfaces in 
Appendix C condition. Later, Anderson and Tsao (Ref. 3) have further supplemented the Appendix C studies 
of Reference 2 with recent data from both SLD and Appendix C tests. 

It was concluded from those two references that acceptable scaling results could be achieved by 
matching the Po, A c , n 0 , and We/,. With scale model size selected, by matching scale and reference values of 
We/ the scale velocity can be determined. By matching p 0 the scale MVD can be found. Reference 2 also 
showed that the effects of temperature and LWC are not independent, but interact through the freezing 
fraction. Therefore, with scale LWC chosen, by matching n 0 the scale temperature can be calculated. Finally, 
by matching A c the scale time can be established. For the scale test, then, only temperature, velocity, MVD 
and time have to be calculated from the known (reference) values of the similarity parameters. 

For swept wing icing scaling, the observations from recent work (Refs. 4 to 7) on ice accretion 
formations on swept wings have suggested that although there are distinct morphological differences in 
ensuing ice shapes on unswept and swept wings, the fundamental physics of ice accretion appears to be the 
same. Thus all the similarity parameters recommended for unswept wing icing scaling should be applicable 
to swept wing icing application provided appropriate modifications to the local collection efficiency and the 
convective heat transfer coefficient in the freezing fraction expression at the stagnation line were made due to 
model sweep angle. 

In a more recent study (Ref. 1) Tsao and Kreeger have shown a best curve fit correlation to calculate the 
Po on swept wings of NACA 0012 airfoil profile. Some experimental evaluations of p 0 were also performed 
for the validation of such correlation. In the present study, a limited number of reference and scale ice shape 
comparisons obtained from NACA 0012 airfoil models at 45° sweep angle were made to assess the proposed 
scaling methods in swept wing icing conditions. 
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Nomenclature 


A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 

c p Specific heat of air, cal/g K 

c p , w .s Specific heat of water at the surface temperature, cal/g K 

d Cylinder radius or twice the leading-edge radius of airfoil, cm 

h c Convective heat-transfer coefficient, cal/sm 2 K 

hj Water film thickness, cm 

ho Gas-phase mass-transfer coefficient, g/sin 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

LWC Cloud liquid-water content, g/m 3 

MVD Water droplet median volume diameter, pm 

n Local freezing fraction, dimensionless 

n 0 Stagnation-point freezing fraction, dimensionless 

p Pressure, Nt/m 2 

p w Vapor pressure of water in atmosphere, Nt/m 2 

p mv Vapor pressure of water at the icing surface, Nt/m 2 

r Recovery factor, dimensionless 

Re fi Reynolds number of water drop, dimensionless 

s Distance along airfoil surface measured from stagnation line, cm 

tf Freezing temperature, °C 

t s Surface temperature, °C 

t Air temperature, °C 

T Absolute air temperature, K 

V Air velocity, kt 

We/. Weber number based on model size and water properties, dimensionless 

Po Stagnation-point collection efficiency, dimensionless 

cp Droplet energy transfer parameter, °C 

X Droplet range, m 

^-stokes Droplet range if Stokes Law applies, m 

A Sweep angle, degrees 

A/ Latent heat of freezing, cal/g 

A,. Latent heat of condensation, cal/g 

p Air viscosity, g/m s 

0 Air energy transfer parameter, °C 

p Air density, g/m 3 

p, Ice density, g/m 3 

p„, Liquid water density, g/m 3 

a Surface tension of water over air, dyne/cm 

x Accretion time, min 
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Subscripts: 

R reference 

S scale 

/ water film 

st static 

tot total 

A at sweep 


Similarity Parameters 

The similarity parameters used in this study were based on the work originally done by Ruff (Ref. 8). 
The scaling method involved matching scale and reference values of the key similarity parameters, p 0 , A c , 
n 0 , and We t . The equations for the similarity parameters will be presented here without much discussion. 
Readers who are interested in the physical descriptions and detailed derivations of these parameters are 
referred to References 2 and 3 and the references given therein. 

To maintain the droplet trajectory similitude, Langmuir and Blodgett (Ref. 9) introduced the modified 
inertia parameter, Ko, defined as 


Kn = 1 

0 8 X 


X 


Stokes 


[k-1 

8 


(1) 


for £>0.125, to describe the inertia of droplets in an air stream flowing around a cylinder of radius d 
positioned normal to the flow. In Equation (1), K is the inertia parameter. 


_ p tl , MVP 2 V 
18 d\i 


Departing slightly from Langmuir and Blodgett in this study, d represents twice the leading-edge 
radius of curvature for airfoils. For the NACA 0012 airfoil model, a leading-edge radius of 0.0158c was 
used (see Abbott and von Doenhoff (Ref. 10)), where c is the airfoil chord. In Equation (1), A./A. Slokes is the 
droplet range parameter, defined as the ratio of actual droplet range to that if Stokes drag law for solid- 
spheres applied. It is a function only of the droplet Reynolds number, Reg. 

FMVDp 

Reg = L (3) 

0 

This study uses a curve fit to Langmuir and Blodgett’s tabulation of the range parameter as given in 
the following expression: 


^Stokes (0.8388 + 0.001483 Re g +0.1847^^) 

of more practical interest than K Q is the collection efficiency at the stagnation point, p 0 , which was shown 
by Langmuir and Blodgett to be a function only of K 0 , 
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(5) 


It should be noted that the (3 0 shown in Equation (5) denotes the collection efficiency at the stagnation 
point for unswept NACA 0012 models. As for a NACA 0012 model at a given sweep angle A (<45°), this 
study uses a curve fit from Reference 1 to calculate the collection efficiency along the stagnation line, p 0 ,A, 

Po,A =Po(cosA) (6) 


Thus the droplet trajectory similarity is satisfied if AT 0 ,s = K 0 ,r and A $ = A« (so is (Po,a)s = (Po.a)«), and the 
scale drop size, i.e., scale MVD, is determined. 

To ensure water-catch similarity, the accumulation parameter is introduced: 


, LWC V x 

A c= 7 

d Pi 


(7) 


If all the water impinging on the leading edge freezes at that location and the leading-edge collection 
efficiency is 100 percent, A c directly becomes a measure of the normalized thickness of ice that will 
accrete. The scale accretion time can be found from A c?s = A c R . When super-cooled water drops strike an 
aircraft surface, they may not freeze immediately on impact. The freezing fraction is the ratio of the 
amount of water that freezes in a specified region on the surface to the total amount of liquid water that 
reaches that region. Thus, local ice thickness depends on Po.A'f and freezing fraction. Because each local 
ice thickness around the model defines the overall shape of the ice, the freezing fraction obviously has a 
major influence on ice shape. The freezing fraction is influenced mainly by the ambient temperature, the 
LWC of the cloud and the aircraft velocity. 

The rate at which the water freezes on a surface depends on the magnitude of local heat transfer 
imbalance. For glaze ice, it is known that the fraction of water that freezes is less than unity, and the 
motion of unfrozen surface water can have an effect on the resulting ice shape. Therefore, it is important to 
maintain surface energy and surface-water dynamics similarities for glaze ice accretions. The freezing 
fraction is formally defined as the ratio of the amount of water that freezes at a given surface location to 
the total amount of water that impinges at that location. From Messinger’s (Ref. 1 1) steady-state surface 
energy balance analysis, the freezing fraction at the stagnation point can be written as 


«n = 


' P,WS 


A 


( 8 ) 


The key terms in this formulation include (|) and 0, which have dimensions of temperature and relate to the 
water drop energy transfer and air energy transfer, and b, the relative heat factor, which was first 
introduced by Tribus, et al. (Ref. 12) 


q> = tf-t s , 


V 2 


(9) 
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Equation (10) from Ruff (Ref. 8) has included compressibility effects. Various incompressible forms 
of 0 have also been used by Charpin and Fasso (Ref. 13) and others; however, the differences are not 
significant mainly due to the fact that, for most icing conditions, the Mach number is relatively low. 

As for the sweep angle effect on the convective heat transfer along the stagnation line of a NACA 
0012 wing section, it was shown by Reshotko and Beckwith (Ref. 14) that for incompressible flow the 
ratio of swept to normal convective heat transfer coefficient became 

h c,A =h c (cosA )2 (12) 


It should be noted that the h c shown above represents the convective heat-transfer coefficient at the 
stagnation point for unswept NACA 0012 models and for brevity its expression can be found in 
Reference 2 for interested readers. Since for most icing conditions the Mach number is relatively low 
Equation (12) is used. Some experimental evaluations of this expression for calculating the freezing 
fraction on a swept wing were made in this study. 

In 1988 Bilanin (Ref. 15) presented a Buckingham- 7t analysis in which he concluded that surface- 
water phenomena had to be included in icing scaling methods. Olsen and Walker (Ref. 16) and Hansman, 
et al. (Refs. 17 to 19) studied surface effects and surface water during ice accretion, presenting additional 
evidence that these were important phenomena to consider in ice accretion. From the close-up photographs 
of these research studies, it was observed that for glaze ice accretion unfrozen water on the ice surface 
tended to coalesce to form beads. These beads sometimes were swept downstream and sometimes froze in 
place. Bilanin (Refs. 1 5 and 20) also argued that drop splashing on impact might affect the shape of the ice 
accreted. 

Hansman and Tumock (Ref. 17) found that when a surfactant was added to the icing spray water, the 
ice shape appearance and shape changed significantly, with the glaze horns moving toward the leading 
edge. Clearly, then, surface tension, and by implication, surface phenomena, have a significant role in the 
physics of ice accretion. 

In 2003 Anderson and Tsao (Ref. 21) had provided experimental evidence from past studies to show 
that a similarity parameter dependent on the ratio V x c'7a : must be included in scaling methodology to 
account for surface-water dynamics effect in glaze ice accretions, although the powers x, y, and z are not 
yet determined. The length may not be chord itself but rather some physical characteristic L related to 
chord; for example, the water-film thickness. Likewise, the velocity could also be of the water-film that is 
related to V. Thus a Weber number based on L and V 


Wei =^P*L (13) 

a 

has been suggested as a potential additional similarity parameter to supplement Ruffs basic scaling 
method. Studies by Bartlett (Refs. 22 and 23) and Oleskiw, et al., (Ref. 24) found no measurable effect of 
pressure on ice shape. These observations suggest that water density is a better choice than air density for 
Equation (13). In this study the We L is based on the twice the nose radius of the airfoil: 
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(14) 


We L = 


V 2 d? w 

a 


with the understanding that Lxd. The scale velocity found from matching We^s = We^j is 

, rt /2 


Vc=V, 


JL 
j s ) 


(15) 


Test Description 

Facility, Model, and Procedures 

The icing tests were performed in the NASA Glenn Icing Research Tunnel (IRT). The IRT is a closed- 
loop, refrigerated, sea level tunnel with a 1.8 by 2.7 m rectangular test section. The icing cloud is 
generated by operating 10 spray bars, a configuration in use since 1998. The IRT cloud calibrations for 
both Appendix C and SLD conditions used for these tests were released in the January of 2010. The LWC 
and MVD measurements were made using methods reported previously (Ref. 25). The LWCs reported in 
this paper are based on an analysis of the LWC calibration data completed in April 2009 for both Standard 
and Mod-1 nozzles. In addition, since only a limited range of MVD-LWC combinations for speeds of 51, 
77, 103, and 129 m/s (100, 150, 200, and 250 kt) have been calibrated to date in the SLD regime, all SLD 
reference cases selected in this study were limited to these calibrated conditions. 

The aluminum test articles used in this study were NACA 0012 Adjustable Sweep Airfoil Models. The 
reference and scale airfoils have a 0.914 m (36 in.) and 0.381 m (15 in.) chord measured normal to the 
leading edge, and a 1.524 m (60 in.) and 0.609 m (24 in.) span, respectively. The reference model was 
mounted in the tunnel on the floor at 45° sweep angle position (left picture in Fig. 1). The scale model was 
mounted in the tunnel on a stand with a small extension that allows pivoting of the airfoil to sweep angle 
of 45° in the mid-span of the test section where the cloud uniformity is the best (right picture in 
Fig. 1). Blue horizontal lines at the leading edge were drawn on the surface of the airfoil at the tunnel 
vertical center to locate ice-tracing templates and to serve as a distance scale in some pictures. Because of 
the quick-start capability of the IRT spray system, no shielding of the models was required during the 
initiation of the spray. 



Figure 1 . — Swept NACA 0012 airfoil models installed in IRT test section. 
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To record the ice shapes, a thin slice was first melted through the ice normal to the model leading 
edge. A cardboard template was then placed into this slit and an outline of the ice shape traced by pencil, 
giving a two-dimensional cross section of the ice. Tracings were taken at the vertical center of the tunnel 
(91 cm from the floor) and at 2.5 cm below the center. The ice shapes so recorded were digitized using an 
automated line-following feature in the image-analysis software, SigmaScan Pro (Ref. 26). The results 
presented in this study are from IRT test entry in April 15-16, 2010. Though the actual swept wing ice 
accretions could exhibit very complex three-dimensional ice shape features like scallops, the shape 
differences between the two tracing locations were never significant. Therefore only centerline shapes will 
be reported here. It should be noted that in general the IRT cloud uniformity, especially for LWC 
distribution, in SLD regime is not as good as in Appendix C for any given air speed. 

Uncertainty Analysis 

Estimates of the uncertainty in the reported average conditions were made by considering inherent 
errors of instruments, temporal fluctuation and spatial variation of the instrument readings in the test 
section, and uncertainty in tunnel calibration of MVD and LWC. Recorded air temperature was believed to 
be accurate to ±0.5 °C, and the uncertainty in air velocity was estimated to be ±1 m/s. For Appendix-C 
conditions the net uncertainty in MVD was estimated at ±12 percent. For SLD conditions it may have been 
as much as ±20 percent. These uncertainties are not referenced to an absolute value of MVD, which is 
unknown. Repeatability and scatter in the LWC calibration data suggests the uncertainty is about 
±12 percent for both Appendix-C and SLD conditions. 

The test-parameter uncertainties were used to estimate the following uncertainties in the similarity 
parameters for the Appendix-C tests: 3 percent in p 0 , 12 percent in A c , 13 percent in n 0 and 2 percent in 
We/,. For the SLD tests the uncertainties were: 4 percent in p 0 , 12 percent in A c , 12 percent in «o and 
4 percent in We/.. 


Test Results 

Due to very limited IRT test time available for the scaling method evaluation, this study only used the 
IRT calibrated SLD drop sizes as reference conditions. These reference cases included MVDs of 44 and 
93 pm, and airspeeds of 100 and 150 knot. A baseline condition of V= 100 knot and MVD = 44 pm was 
chosen, and air temperature was varied to create a range of stagnation-point freezing fractions (n 0 ). Scale 
conditions were then calculated by matching p 0 , «o. A c , and We L with selected scale LWC value for each of 
the reference case. For each figure, reference ice shape was shown shaded, while a solid line indicated the 
scale ice shape. The table below each figure gave the test conditions and similarity parameters for each 
pair of reference and scale tests. The conditions given were the average conditions recorded over the 
duration of each test, which can sometimes differ slightly from the planned set points. The parameters in 
the tables were calculated from these average conditions. 

V = lOOkt (16) 

For the baseline condition of MVD = 44 pm, there were three sets of reference and scale ice shape 
comparison for scaling method evaluation. 

Stagnation-point freezing fraction of 0.41 . — The first pair of ice shape comparison was shown in 
Figure 2, and it was for stagnation-point freezing fraction of 0.41 . The scale ice shape did simulate the 
reference main ice shape and feather region reasonably well, even including smaller feathers farther aft on 
the surface. Further examining the front, side and close-up view of reference and scale ice shapes (Fig. 3), 
it was noted that the scale ice scallop tip structure was somewhat less pronounced than the reference ice. 
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Stagnation-point freezing fraction of 0.57. — The second pair of ice shape comparison was made for 
stagnation-point freezing fraction of 0.57. In Figure 4, it was shown that the scale ice shape was able to 
duplicate the reference main ice shape and feather region quite well but the size of reference ice was 
clearly off and under-sized. This was believed to occur because the slits were not totally clear from ice and 
the cardboard template was not able to go all the way in as freely as it should be. Consequently the 
resulting ice shape tracing became undersized. This mishap could happen rather often if one did not take 
enough time to melt the ice completely prior to the tracing. Unfortunately there was not any simple means 
to recognize this problem until the ice shapes were digitized and compared. Further examining the front, 
side and close-up view of reference and scale ice shapes shown in Figure 5, it was noted that the scale ice 
scallop tip structure also simulated the reference ice scallop tip structure fairly well. 

Stagnation-point freezing fraction of 0. 74. — The third pair of ice shape comparison was shown in 
Figure 6 for stagnation-point freezing fraction of 0.74. The scale ice shape duplicated the reference main ice 
shape and feather region very well, even including smaller feathers further aft on the surface. Further 
examining the front, side and close-up view of reference and scale ice shapes shown in Figure 7, it was noted 
that the scale ice scallop tip structure also simulated the reference ice scallop tip structure fairly well. 

Stagnation-point freezing fraction of 1.0. — Figure 8 showed the last pair of ice shape comparison for 
baseline condition at stagnation-point freezing fraction of 1 .0, i.e., rime ice shapes. It was shown that the 
scale ice shape was able to simulate the overall reference rime ice shape features which included the main 
ice shape and feather region quite well except the scale ice shape size was larger than the reference ice 
because a similar tracing template problem shown earlier in n Q = 0.57 case has occurred again. Typically if 
the cloud’s LWC were off from what it should be from the scaling calculation, the resulting ice shape not 
only would have different size but also would have quite different ice shape features. The size mismatch 
observed in the figure was clearly not of this nature but an operational mishap. Further examining the 
front, side and close-up view of reference and scale ice shapes shown in Figure 9, it was noted that the 
scale rime ice structure also simulated the reference rime ice structure fairly well. 

V = 150kt (17) 

For this higher reference velocity, there were two sets of reference and scale ice shape comparison 
obtained for scaling method evaluation. One was for the baseline condition of MVD = 44 pm and the other 
was for larger SLD drop, i.e., MVD = 93 pm. 

Drop MVD of 44 pm, stagnation-point freezing fraction of 0.5. — Figure 10 showed the reference and 
scale ice shape comparison for constant We t method. The scale ice shape simulated the reference main ice 
shape and feather region well. It was noticed during the test that some large feathers neighboring the 
reference main ice were accidently melted away while making the tracing. As a result it made the reference 
ice shape appear to be a bit smaller than the scale ice shape. If one examined the front, side and close-up 
view of reference and scale ice shapes as shown in Figure 1 1, it was shown that the scale ice’s complete 
scallop tip structure did duplicate quite well the reference ice complete scallop tip structure. 

Drop MVD of 93 pm, stagnation-point freezing fraction of 0.51. — Figure 12 showed the reference and 
scale ice shape comparison for this large drop size of 93 pm at freezing fraction of 0.51. The scale ice 
shape also simulated the reference main ice shape and feather region well. If one further examined the 
front, side and close-up view of reference and scale ice shapes as shown in Figure 13, it was shown that the 
scale ice’s complete scallop tip structure did duplicate quite well the reference ice complete scallop tip 
structure. 
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Figure 2. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 100 kt, MVDr = 44 pm; no = 0.41. 



Figure 3. — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 2. 
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Figure 4. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 100 kt, MVDr = 44 pm; no = 0.57. 



Figure 5. — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 4. 
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Figure 6. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 100 kt, MVDr = 44 pm; no = 0.74. 



Figure 7. — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 6. 
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Figure 8. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 100 kt, MVDr = 44 pm; no = 1.0. 



Figure 9. — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 8. 
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Figure 10. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 150 kt, MVDr = 44 pm; no = 0.5. 



Figure 1 1 . — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 10. 
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Figure 12. — Reference (shaded) and scale (solid line) ice shape comparison for Vr = 150 kt, MVDr = 93 pm; no = 0.51 . 



Figure 13. — Front, side and close-up views of reference (U) and scale (L) ice 
shapes of Figure 12. 
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Summary/Conclusions 

Limited icing scaling tests were performed in the NASA Glenn Icing Research Tunnel to assess the 
modified scaling methods for swept wing icing application. The reference tests used a finite span 
9 1 ,4-cm-chord adjustable sweep airfoil model ofNACA 0012 profile at velocities of 100 and 150 knot and 
MVD of 44 and 93 pm. Scale-to-reference model size ratio was 1 :2.4. All tests were conducted at 0° angle 
of attack (AoA) and 45° sweep angle. 

The preliminary result showed that good agreement of both size and shape of the main ice accretion 
and feather regions was observed. The present SLD reference tests were made with velocities of 100 and 
1 50 knot, and these observations may not be valid for higher velocities. Additional test is required to 
identify the range of test conditions, both Appendix C and SLD regime, in which the modified scaling 
methods for swept wing icing may not apply. 

Also it was evident from the ice shape comparison presented in this study that the traditional ice 
tracing method to document straight wing ice shapes was not sufficient to properly characterize the highly 
three dimensional ice shape features, e.g., the complete scallop tip, that are prevalent in swept wing ice 
accretion. Currently some effort is taken by NASA Icing Branch to develop a three-dimensional ice shape 
characterization method for swept wing icing. 
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